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■ Abstract 

■ We analyzed T-violation in neutrino oscillation by using perturbation methods with re- 
spect to Am2iL/2E and 6a{x)L/2E, where 6a{x) represents the matter density fluctua- 
tion from its average value. We found that the matter contribution to T-violation arises 
from interferences between Aml'^^L/2E and 6a{x)L/2E. In the 2nd order, the symmetric 
and asymmetric matter density fluctuations give effects to the sin 5 (intrinsic) and the 
cos 5 (fake) parts of T-violation. We give their analytic forms and analyze the matter 
contribution to the sin 5 and cos 5 terms. We found that, for L = 3000km, both the 
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symmetric and asymmetric matter density fluctuations give negligible contributions to 
T-violation, and that thus the constant (average) matter density gives a good approxi- 
mation. On the other hand, we argue that, for L = 7000km or longer length, T-violation 
turns out to become very small due to cancellation between the 1st and the 2nd order 
terms. This shows that the constant (average) matter approximation is not valid. 

1 Introduction 

A high-intensity neutrino source based on a muon storage ring, which is now generally 
called a neutrino factory[l], has attracted growing interest from theorists and experimen- 
talists[2]. One of the important physics potentials at neutrino factories is to measure 
a possible non-zero CP violation phase {S) in the 3-generation neutrino mixing matrix, 
the Maki-Nakagawa-Sakata (MNS) matrix[3]. To measure the phase 6 in the neutrino 
oscillation, one way is to compare the CP-conjugate oscillation processes, P(i^a — J^/s) 
and P{Va — * V/s), and the other way is to compare the T-rcvcrscd oscillation processes, 
P{^a ^p) and Pii'p — > Va)- The study of CP violation has been extensively studied[4]. 

T-violation has not been discussed seriously so far. It is mostly due to the difficulty 
of identification of v^. appearance at a neutrino factory, since the detection of wrong- 
signed electrons is hard. However, there have been some attempts to disentangle wrong- 
signed electrons which enables us to search for T-violation at a neutrino factory[l]. Its 
experimental feasibility studies are now undertaken. One of the advantages of the search 
for T-violation is to expect relatively small contribution from matter, whereas in the 
search for CP violation the fake CP-odd effects from matter dominates over the intrinsic 
CP violation for a long baseline length (such as more than a few thousand km) and 
thereby the measurement of the intrinsic CP violation becomes challenging. 

T-violation in matter arises from the intrinsic contribution with matter modification, 
which is proportional to the CP phase, and the fake matter contribution. The constant 
matter density gives an effect to the intrinsic T-violation which is proportional to the 
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sin 5 term, which has been discussed extensively [5, 6, 7]. However, up to now, it has not 
been addressed about possible contributions from the symmetric and also asymmetric 
matter fluctuations deviated from the average constant density. 

In this paper, we consider the symmetric and asymmetric matter fluctuations from 
the average density and treat them in the perturbation method developed by Koike and 
Sato[8], and Ota and Sato[9], where the quantities, AmliL/2E and Sa{x)L/2E, are 
considered as perturbative Hamiltonians, which may be small for most of the cases. The 
average matter density is included in the unperturbed Hamiltonian, so that the constant 
matter contribution (the average matter) was taken into account. The matter fluctuation 
is separated into the symmetric and asymmetric terms, Sa{x)sL/2E and 5a{x)aL/2E, 
respectively. Ota and Sato considered the symmetric part by using the preliminary 
reference earth model (PREM)[10] and analyzed the 1st order corrections for Aml-^^L/2E 
and 6a{x)sL/2E. In the 1st order, T-violation arises only through Am2iL/2E. The 
matter fluctuation does not give any contribution to T-violation in this order. 

We examined T-violation in the 2nd order perturbation of Am|iL/2£' and 5a{x)L/2E. 
Our motivation is to obtain the analytic expression of T-violation in the 2nd order and 
to examine how large the 2nd order contribution from the symmetric and asymmetric 
fluctuations is. We found that the interference term between Am2iL/2E and 6a{x)L/2E 
gives some contributions to T-violation. In particular, the symmetric matter fluctuation 
contributes to the sin 5 part, while the asymmetric matter fluctuation does to the cos 5 
part. We estimated these contributions for L — 3000km case and found that these 
contributions are negligibly small for most energies. As a result, the constant matter 
approximation works well. On the other hand, for L = 7000km or longer, the 2nd order 
term becomes comparable or larger than the 1st order term, so that the matter fluctuation 
can not be neglected and the constant matter density approximation fails. 

In Sec. 2, the perturbation formula is given, and the 0th and the 1st order contributions 
with respect to /\m^^L/2E arc given in Sec. 3. The general discussion on the contribution 
from the asymmetric matter profile to T-violation is given in Sec. 4. In Sec. 5, the 2nd 
order contribution from Am\iL/2E is given. The interference term between A.m\]^L/2E 
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and 5a{x)L/2E is presented in Sec. 6 by assuming the linear dependence for Sa{x) and 
the numerical analysis of these interference terms is given in Sec. 7. The summary is given 
in Sec. 8. 



2 The perturbation formula 

The formula to evaluate the neutrino transition probabilities perturbatively with respect 

to the small quantities, Arn^iL/lE and 6a{x)L/2E, has been developed by Koike and 
Sato [8], and Ota and Sato [9]. Ota and Sato used this formula to estimated the 1st 
order terms of Aml-^L/2E and the symmetric matter fluctuation, 5a{x)sL/2E. Here, we 
calculate the higher order terms with respect to the symmetric, Sa{x)sL/2E, and the 
asymmetric terms, 5a{x)aL/2E. Firstly, we outline their method. 
We begin with defining the neutrino mixing matrix as 

U = e^^^^Miag(l,l,e^^)e^^^^5e^^^^2 



SxCy CxSyS z^- CxCy ^x^y^ Sy^z^ 



i5 

\ ^x^y ^x^y^z^ ^x^y ^x^y^z^ ^y^z^ J 

where \j {j — 2,5,7) are Gell-Mann matrices and = cos^^ and = sin^^. The 
angles 9x, 9y and 9z correspond to 9i2, ^23 and ^13, respectively, where % are defined 
in the particle data group [11]. Since the Majorana CP- violation phases are irrelevant 
to the neutrino oscillations (flavor oscillations) [12], we neglected them. If we multiply 
the irrelevant phase matrix diag(l, 1, e~*'^) from the right-hand side of U, we obtain the 
standard form[ll]. The relation between the flavor eigenstates, 1^^) (a = e,ii,T), and 
the mass eigenstates, ji^j) {i — 1, 2, 3), is given by 

Wa) = Uaill^i) . (2) 

The evolution of the flavor eigenstates in matter with energy E is given by 

^^|i^^(a;)) = H{x)f3aMx)) , (3) 



(1) 



where Hamiltonian H{x)pa is given by 
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Here Amfj = ml —rn?j with being the mass of |z/j), Gj? is the Fermi couphng constant 
and 

a(.) . 2V2G.,U.)S ^ 7.56 x 10- (^) (^) (^) eV^ , (5) 

where ne(x), Yg and p{x) are the electron number density, the electron fraction and the 
matter density, respectively. For the electron fraction, we use Y^. = 0.5. 
We separate the matter density fluctuation from its average a, 

5a{x) = a{x) — a , (6) 

and consider the deviation 5a{x) as a perturbative term. That is, we solve the evolution 
equation by treating Sa{x)L/2E and AmliL/2E as perturbative terms, because they are 
small for most of the cases of planned neutrino factories. The validity of this perturbation 
was discussed by Ota and Sato [9]. They in fact showed that the transition probability 
of the neutrino oscillation is well reproduced if the 1st order perturbation with respect 
to the symmetric matter profile, 5a{x)sL/2E, is taken into account, where it is assumed 
that the symmetric matter profile is well approximated by the preliminary reference earth 
model (PREM)[10], for L = 3000km, L = 7332km and L = 12000km. 

(a) The definition of Hamiltonian 

Following the work by Ota and Sato [9], we divide H{x) into the unperturbed part 
Hqo and perturbed parts, ifoi and Hi 

H^Hoo + Hoi + Hi{x) , (7) 

where 

Hoo = 2^e^'^^Miag(l,l,e'^) 
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For the later calculations, we use the following quantities; 



/ C5 S5 \ 



tan 26*2 



C/o = e^^«^Miag(l,l,e'^)e^^^-^« 
S2z{Aml^ - Aml^sl) 



SyS^G 



V C^y^z^ ^y ^y^z^ j 



We also define 



With Oi (i = ±,0), we define 

kx — ao — a- , k2 — a+ — qq , 
k — A;2 + /ci = a_|_ — a_ , 



and 



l^-ia.x Q Q >^ 

e-^"o^ 
V e-^"+^/ 



(8) 



C2^(Amii - Amiis^) - a ' 

- (Am^, + Am\^sl + a 

±V{c2.(Amii - Amiis2) - a}^ + sUAml, - Am\s^:f^ . (9) 



(10) 



(11) 



(12) 
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(b) Interaction representation 

In the interaction representation, the interaction Hamiltonians are given by 

Hoi{x)i = e'^°°''Hoi{x)e-'"°''' , Hi{x)i = e^^°°^ii'i(x)e-^^°°^ , (13) 

and the wave functions are presented by = {e'^""'')ai3Wi3ix)) with |fa(0))7 = 

|fa(0)). The evolution equation of li'ai^))! becomes 



The solution is given by 

\i^a{L))i = T exp 



[-i) dxHoi{x)i + Hi{x)i 
Jo 



(14) 



(15) 



a/3 



where T means the time ordered product. Then 



KiL)) 



-iHooL 



-iHooL 



al3 
L 



ML))i 



l + I dxHoi{x)i+ {-i) I dxHi{x)i 
Jo Jo 



+ dx dyHQi{x)iHoi{y)i+ {-if dx dyHi{x)iHi{y)i 

Jo Jo Jo Jo 



+ (-i)2 / dx r dy{Hoi{x)iH^{y)i + H^{x)iHoi{y)i) + ■ ■ ■ 

Jo Jo n 

J 7/3 

= (-S-oo + -Soi + 5f ) + ^-oLoi + -Sf ) + 5oi,i + • • k/3(0)) . (16) 
The transition probability from one flavor eigenstate a to another (5 is given by 

— |('S'oo)a/3r 

+ 2Re (>S'oo)q,(3(5'oi)*^ + (5'oo)a/3(5'i '')*^ 
+ |(^oi)a/3r + \{S'^x\p? + 2Re \{S^x\AsV)lt 



+ 2Re 
+ ••• . 

T- violation is defined by 



(>S'oo)a/3(>S'oi,Ol)a/3 + {Soo)ali{S\ ^fafj + (>S'oo)a/3 (>S'oi,l)a^ 



[17) 



(18) 



and is evaluated with use of the probabilities defined in Eq.(17). 



3 The 0th and 1st order contribution from Hqi{x) 

Ota and Sato[9] calculated the 1st order contribution from Hqi{x). Here, we give a brief 
derivation of their results which are needed to discuss the higher order calculation, 
(a) The S-matrix (0th and 1st order) 

The unperturbed Hamiltonian Hqq is diagonalized explicitly and is given by 



i^oo = t/odiag(a_, Oq, a+)Ul , 



(19) 



where a± and Gq are defined in Eq.(lO). 

The S-matrix for Hqq is easily obtained as 



^00 = e-'^°°^ = UoP{L)Ul 



/ 0+ - C25( 



SyS2z^ 



S2y {cjh- - Si01-) 



CyS2ze 

S2y (ci02- 



1-)/ 

(20) 



where 0± and 0j± {i = 1, 2) are defined in Eq.(ll). The 1st order term of Hqi defined in 
Eq.(8) is given by 

/ Ai \ 



where 



Soi = Uo Aoi Boi 

V Boi J 

( CyVox 

Cj/'Poi —S2yC8Q-m 



i&c2 „i(5^2 



4)Q, 



01 



01 



S2yC5Q^ 



01 



(21) 



7^01 
Ai 



CzAoi + s^Boi , Qoi 

52a;Cz-2 , 



SzAqi - CzBqi 

i3, - —''i 



^01 



4£; ko 



(22) 



(b) The oscillation probabilities (0th and 1st order) 
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The oscillation probability in the 0th order P^°^\ua — > i^/j) is given by |(-S'oo)a/3p, 



2 



'2i sm — , 



= sly Usu.' ^ + clsin^ _ ,|elsin^ ^| , (23) 



with A;, /ci and defined in Eq.(ll). 

The probability for antineutrinos P*^°°)(i^ — > i^g) is obtained by taking 5^—5 and 
a — > —a. 

For the higher order terms, we consider only the contribution to T- violation. The 1st 
order term from Hqi, i.e., the AmliL/2E term is solely by the CP violation phase S and 
is given by 



( A nT \ ^''''21 



sm sm sm — . (24) 

2 2 2 ^ ^ 



where we used sin /ciL + sin A;2-L — sin /cL = 4sin(A;iL/2) sin(A;2-^/2) sin(A;L/2), which is 
proved by using k — k\^ki. This formula includes the constant matter effect. 
Similarly, we find 



which is valid for the constant media as stated in the paper by Krastev and Petcov[5] . 



4 The contribution from the matter term, B.\(x) 

In this section, we give the general formula to evaluate the n-th order effects of matter, 
i.e., the {5a{x)L/2EY order. Then, we evaluate the contributions of up to the 3rd order 
and discuss the general properties of the effects. 

(a) S- matrix elements 
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The interaction Hamiltonian with matter in the interaction representation is expressed 



by 



_ 6a{x) ~ 



f cj sscse 


2 





Then, the n-th order matter perturbation is given by 
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Therefore, the S-matrix is written in general by 
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AS o(-) o2 (-) 

\e Cy/Sijj^ ^y^y^l,n ^y^l,n I 



(28) 



where 



= + ^{"^±(c2.-(£:?^-^J'^^)+s2,(ci"^ + pi"^)) , 



-in) 



(n) 



,(n) 'n(n)\ 



(29) 



Below, we show their explicit forms, which are needed to evaluate the matter effect 
to T-violation up to the 3rd order. The 1st order terms are 



CI 



(1) 



(1) 



^ ^ Jo 2E 

^ 6a{x] 



-ikx 



Jo 



dx- 



£(1) = (-i)e-'"-^ / 



2E 

L 



, 8aix) 9 

fir - -r- 

(1) _ l_,\^-ia^L r ^J<^) 



dx- 



2E • 



(30) 
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The 2nd order terms are 



6? = (-,)^e-"-jf'd./;*^^M^.,c,(.3e-- + c3e-'»), 

ri'' = /; dy'-^'-^sl (4 + 4.'"-') . (31) 

Finally the 3rd order terms are 

^ ^ ^ Jo Jo Jo 2E 2E 2E 

^ ^ ' Jo Jo ^ Jo 2E 2E 2E 

xsscs (cle^^^^ + slc1{e'''y + e^'=(^-2'+^)) + sje''"') . (32) 

It should be noted that Vi^^ and J-'i"'^ are derived from Ci"-* and S["'\ by exchanging 
between Sz and c^, and also a_ and a+. 

In order to examine the general properties of these quantities, we divide the matter 
fluctuation Sa{x) into the symmetric part, Sa{x)s, and the asymmetric part, 5a{x)a, and 
then expand them in terms of Fourier cosine series, 

'TTT) 

daix)a = E«2n+ie-*^^"+^% gn = — , (33) 

n ^ 

where a„ is a real number satisfying a_„ = a„. Here, we excluded n = in the expression 
of Sa{x)s, because 6a{x) is the deviation from the average of a{x). Below, we consider 
the both profiles together. 
The 1st order terms are 

rW , -r)(l) _ «25 / a2n 

£;(!) ^ JF« = 0. (34) 
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The 2nd order terms are 



"S22C22 



E 



^(2) _ ^|2) 



"S22C22 



E 



(2^)' L+t^odd + ^n){k + gn + g„0 



(t>+ , 



"'22 



8(2^)^ 



f /y^ ((-ir + l)((-ir + l)ana^ \ 



(35) 



Finally, the 3rd order terms are 



2(2^)5 



^2 -l)anan^a^ 

^^"^ + ^")(^" + 9" + ^rn){k + qn + qm + Ql) 



s2 



n,m,l 



{k + n {k + Qn 



(k + qi) ^_^{k + qn){k + qm) 



4 ^ 



E 



{k + qi) n {k + qn 



i(f)- 



(36) 



(b) The general properties 

We first discuss the n-th order contributions of the matter to the transition probabil- 
ities, which are obtained by computing 

P^''''\i^a^r^p) = ReUsoo)ap{Si''^*)ap+ E isi%fs{st^*) \ . (37) 

\ l+m=n 



Since the (-Sj"^)^^ = {SD we conclude that 



(38) 



Therefore, the matter itself does not give any effect to T-violation for 
For other channels. 



Vt channel. 



(^e-^.) = ^Re I 252,0-/3} J* + E f^uA 

l+m=n 



(+)^(+)* 
m 



(39) 
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and P(^'"'^(f^ — > Ue) is obtained by changing the superscript (+) to (— ). The transition 
probabihties for u,. Vr and Vr ^ are obtained from those for i/^ — > i^ij, and — > fg, 
by changing the coefficient with c^. 

We consider the 1st order of the asymmetric matter fluctuation, which contributes 
only C'l^ — V^l\ in contrast to the C'l^ + V^^^ which is due to the symmetric matter 
profile. In addition, Si'^ — and Ti'^ — 0. Therefore, the asymmetric matter fluctuation 
gives a\^i = 0, which shows that it does not contribute to Vt channel. Also we find 
/Sj^'* = ±(Ci^'* — I^i^''), which is proportional to 0+. Since 

= -2i sin kL , (40) 

which implies the vanishing contribution. Thus we conclude that the asymmetric matter 
profile does not contribute to the transition probability in the 1st order. It contributes 
to the transition probability in the 2nd order. 

(c) T- violation only from the matter 

(c-1) The 1st order effect 

The 1st order contribution to AP^^^ is from the interference between 5*00 and S^^ as 
we see in Eq.(17). In this case, the C^^ — V^^^ in S^^ contributes to T-violation. Since 

—T)^x^ is proportional to 0+, while 5*00 to T-violation is proportional to Re[0_0y , 
which is zero. Therefore, there is no 1st order effect to T-violation. 

(c-2) The 2nd order effect 

The 2nd order effect comes from and the interference between 5*00 and S^^ . 

Similarly to the 1st order contribution, Re[0_(Ci^''* — V>y^*)\ in Sf^ contributes to the 

(2) (2) 

asymmetry. Since C\ —T)\ is proportional to this term does not contribute. As 
for l^f term, Re[(£f ^ - J^J'^)(Cf ^* - ©i'^*)] and also Re[(Cf ^ + PS'^)(Cf ^* - ^*)] 
contribute. Since C^^ ^v'^^^ is proportional to 0_ and —T'^^ consists of 0_ and 
they do not contribute. As a result, there is no 2nd order effect. 

(c-3) The 3rd order effect 

The 3rd order effect comes from and S^^ S^^* . Since — "D^^ contains 
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0+ and i(f)-, no contribution arises from 5'oo'S'i term. Similarly, for terms, the 

terms contributing to T-violation are Re[{S[^^ - J^f ^) (C{^^* -pf^*)] , Re[(Cj^^ +V^^^){Ci^^* - 
T>^^*)] and terms that the superscript 1 and 2 are exchanged. Again it is clear that these 
terms do not contribute. 

We conclude that there is no contribution to T-violation through the matter terms 
up to the 3rd order, i.e., {6a{x)L/2E)^ terms. The vanishing of the 3rd order term is 
sufficient enough to deal with the actual situation. Thus we do not pursue the further 
investigation, although we expect that the higher terms will not contribute either. 



5 The second order contribution from Hi 
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We consider the (Am|iL/2£^)^ effect. For this, there are two contributions. One is from 

\{Soi)pa\'^ and the other is from 2Re[(S'oo)/3a(5'oi,oi)^c«]- 
We first compute 5'oi,oi defined in Eq.(17) 



Jo Jo ^ 

The quantity of the bracket becomes (Am|]^S2a;/4-E)^ multiplied by 

^ ^2 _^-iki{x-y) 



Then, we find 



„2 „iki{x-y) 



+ K--z^' 



-ik2{x-y) 





g2 _^ik2{x-y) ^ 



01,01 



where 



01,01 







Coi,oi \ 



/01,01 







\^^01,01 

/ "01,01 

J8^ r{-) 



^ 31,01 / 

"'?;Poi,01 



e*'^Sy/?01,01 SyOioi.Ol + 2c^^oi,01 SyCyftoi/oi " S2yQo\,^ 



e *'^Cj^/3oi,oi 
(-) 



01 



V e"^Cj,/3oi,oi ^yf^yf^WfM ~ ^2yQo\,o\ c^Q;oi,oi + 2s^^oi,oi / 



„(-) 



,2^,{-) 



(41) 



(42) 



(43) 



"01,01 
Poi,oi 



(^01,01 + -^01,01) ± (C25(^01,01 — -^01,01) + 522(^01,01 + 1^01,01)) , 

—•§22(^01,01 — -^01,01) + C22(Coi,oi + "^^01,01) i (Coi,oi — "^^01,01) ■ 



(44) 
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We compute the 2nd order (Am|iL/2£')^ contribution to T-violation. Since (5'ol,ol)/^T = 



(5, 



01,01 Jr^, 



there is no contribution to T-violation for the v,, and Ut oscillation channels. 



For the and i/^ channels, the Coi,oi — ^^oi,oi contributes because 

2Re[(5'oo)/ie('S'oi,Ol)^e ~ (>S'oo)e/i('S'oi,Ol)e/i] = "SyS2iRe[0_ (Cg^ — I^oi,Ol)] • 

The same holds for the and Vr channel. 

Coi,oi and ©01,01 are given by 

<^oi,oi — 



(45) 



•^01,01 



;S2{z-z) 



:S2(z-z) 



1 



kki 



+ 



1 



kik2 



(46) 



Since 



01,01 



'Aml^S2xY 1 



01,01 



-S2(z 



"21-: 

V 4£; 

/ Aml^S2xY 1 



kk2 



02- 



1 1 

+ 



.fcA;2 fciA;2 ^^1 

because 02- + 0i- = <P- and k2 — k — ki, there is no contribution to these channels 
The |(5'oi)/3aP also give a null contribution because 



(47) 





\{Soi)fMef - 


\Cy{CzAoi + S^fioi)!^ - \Cy{CzAoi + S^Bqi) 




|('S'oi)err ~ 


('S'oi)rel^ — 


\sy{c^Aoi + s^Boi)!^ - \sy{c^Aoi + s^Boi) 








\{e''sl-e~''cl){s,Aoi-c,Boi)\' 








-\{e-^Ssl-e^'4){s,Aoi-c,Bo,)\'^0. 


(48) 



Therefore, we find no effect from the 2nd order term of Hqi, 



A pr(01,01) ^ A pT(01,01) ^ A pT(01,01) ^ Q 

t^et^^i i^el^T t'/ii^T ■ 



(49) 



6 The interference between Hqi and Hi 

In this section, we consider the {Am\^L/2E){5a{x)L/2E) contribution to T-violation. 
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The S-matrix is given by 





( ° 




] 






>^'oi,i 











I 


Boi,l 


j 





(50) 



where 



^'01,1 — 

Boi,i = 



AE 
^ AE 



" 

(5a(?/) 



ik\x 



2£; 

[ dx [ dv ^'^^^'^ Sz (sz 
Jo Jo 2E ^ 



4S 



:e-^"°^ / dx r dy^-^s-z (-s~zSz-ze-''^^ + 
Jo Jo 2E ^ 



Z'-'Z — Z^ ' ^Z^Z — Z^ j 1 

—i{ky—kix)^ 



Then, the S-matrix is given with use of v4.oi,i and i3oi,i by 

/ CyVoi,l -SyVoi,! \ 



'S'01,1 



where 



(51) 



, (52) 



7^01,1 = (civ4.oi,i + SzBoi^i) , Qoi,i = (siAi.i - CzBoi,i] 

^01,1 ~ ('^2Al,l + -^z^Oljl) ' 2oi,l ~ ('^zAl.l ~ '^zBqi ij 



(53) 



(a) The contribution from the symmetric matter fluctuation 

The interference terms between Am|^L/2£' and the symmetric matter fluctuation, 
Sa{x)s, contributes only for the term proportional to the CP phase, sin 5. There are 
two contributions. One is from the ^'oi^'i^''* and the other is S'oo^S'oi ig terms. We take 
the PREM[10] as a symmetric matter. Ota and Sato expanded the PREM in the cosine 
series[9], which is expressed as 5a{x)s in Eq.(33) with appropriate values of coefficients. 
The Si^J* is given by Eqs.(28) and (34), where we take only the even n case, i.e., from 
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C[^^ +V[^\ We find 



E 



+ 



C2z 



E 



(l2nk 



. kiL . /C2I/ . kL 



„ , I sm sm sm 

.2 _ ^2 ^ I 2 2 2 



The term 5*01, is is obtained once we compute Aqi^u etc., which are given by 

^fn2iS2x f f CzCz—z ^z^z 



«4oi,ls — Aq 



01,1s 



4E 



0'2nkl 



01,1s 



B' 



/ci k2 



k 



E 



a2nk 



2 n=l,2,. 



01,1s 



X 



0'2nk2 



02- - 



E 



^1 \n=ti.- E{k' - qL) 



0- 



Now, we find 



2Re[(5oo)eM(^oi,is):^ - (e ^ //)] 

_ ^^21*2x^22/52555 / Cz-z^z Sz-zS'z 



X 



E 



0'2n ( c\ki 



ki 



+ 



k,. 



S'k2 



kiL . k2L . kL 



n=l,2,-- V^l y2n ^^2 y2n 

Thus, we obtain the sin b part of T- violation as 



2_.2 ^2_„2 ;sm — sm — sm-^ 



'-S2xS2yS2zS5 



^zCz—z ^ ^z^z—z 



ki 



ko 



sm sm sm 



x<l+ J2 

n=l,2,- V^l ~ Q2n ^2 ~ ?2n ~ ?2n. 



2 _^ C2i/C 



Similarly, we confirmed that the relation 

holds in the order of {AmliL/2E){6a{x)sL/2E). 
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(b) The contribution from the asymmetric matter fluctuation 

The interference term between Am2iL/2E and the asymmetric matter fluctuation, 
6a{x)a gives the cos 5 part of T- violation. In order to estimate this contribution, we 
consider an extreme case, the linear fluctuation, as a typical asymmetric matter profile, 

6a{x)a = Ota ^ — TT^ ' ^^^^ 

where a represents the fraction of the asymmetric matter profile. 

Firstly, we compute SqiSi^* term, where S'l}* is the 1st order term from the asym- 
metric matter profile. For this, we estimate C'l^ — V^^^ defined in Eq.(30) from 5a{x)a- 
We find 

Now the contribution from the interference term between S'oi and S^^ is obtained by 
using Eqs.(21) and (28) as 

2Re[(5oi)e^(^i):^ - (e ^ //)] 

_ aa /^mlyS2xS2yS2zC& f Cz-zCz Sz-zSi 

2EL E \ ki k2 

(I L kL\ . kiL . k2L . kL 
X — cot — sm — — sm — — sm — . (61) 
\k^ 2k 2 J 2 2 2 ^ ^ 

Next we consider the <S'oo'S'oi,ia contribution. By performing the integrations, we find 
aa AmliS2x ( { CiC^^z . SiSz-z\ ( L 



4EL AE \\ h k2 J V^i H 



H ; 1-01+ - ^7^0l 



Sz^z-z f L .2 

^01, la = — •6'oi,la 



aa Am|iS2a: / ( Cz^z-z , SiS^_i\ I L 



AEL AE 



^z^z—z f L ^ .2 



ki \k 
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Now T- violation for the Ve and oscillation channels is given by 

2Re[(5oo)e/.(5oi,ia):^ - (e ^ //)] 
_ aa Aml]^S2xS2yS2zC5 (Cz-zCz Sz-zSz\ . kiL . k2L . kL 



H ; sm sm sm ■ 



2EL E \ ki k2 J 2 2 2 

f f4 sl\ Lf4 kiL si k2L\\ 



By adding these contributions, we obtain the asymmetric matter effect to T-violation, 

aa Aml.S2xS2yS2zC5 ( CzCz-z , SzSz-z\ . kiL . k2L . kL 



sm sm sm 



cs 2EL E \ ki fca / 2 2 2 

\\k^ kl kl) 2\k 2 ki 2 k2 2 J j 

(64) 

Similarly, we confirmed 

From Eqs.(58) and (65), we find that T-violation is independent of flavor in the 2nd 
order perturbation. 



7 Numerical analysis 

In order to analyze the matter effect, we use the PREM profile for the symmetric matter 
density distribution. 

(a) L = 3000km 

(a-1) The T-violation (intrinsic) asymmetry 

Firstly, we estimate how large the T-violation asymmetry is. In Fig.l, we plotted 
{APj^^^)sj2P{i'e I'fj) as a function of the neutrino energy E, with S — 7r/4, S2x — 
S2y = 1, = 0.1. Here, we used our formula in Eq.(57) for (APj^^^)^^ and the 0th order 
term for P{ve v^). For the symmetric matter profile, we use p = 3.31602g/cm^ and 
the Fourier coefficients, 02, 04, oe and as are —0.045, —0.048, —0.047 and — 0.044g/cm^ 
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respectively, which are derived by Ota and Sato[9] from PREM. We neglect the n > 10 
terms. As we can see from this figure, we can expect about 4 ~ 10% effect for E > 5GeV. 

(a-2) The matter-modified T-violation in the symmetric matter profile 

In the 2nd order perturbation, the symmetric matter gives the contribution to the 

sin 5 of T-violation and the combined formula is given in Eq.(57). 

In Fig. 2, the comparison between our result and the vacuum case of T-violation, 

(^^-^-J vacuum = -^2.S2yS,,C,Ss Siu siu siu . (66) 

is made, for (a) 3GeV< E <30GeV and (b) lGeV< E <3GeV. The values of parameters 
are the same as those in Fig.l. The sohd hue shows the vacuum case and the dash- 
dotted line shows the matter-modified T-violation with the addition of the 1st and the 
2nd order terms. The 2nd order contribution is shown by the dotted hue, but it is hard 
to see because it is almost zero in this scale. That is, the 2nd order term is negligibly 
small and we can safely use the formula given by taking the average density. There is the 
matter enhancement around E — 6GeV which is consistent with the discussion by Parke 
and Weiler[7]. 

(a-3) The matter-modified T-violation in the asymmetric matter profile 

In Fig. 3, we plotted (AP)c^, the fake contribution to T-violation from the asymmetric 
matter fluctuation, for (a) 3GeV< E <30GeV and (b) lGeV< E <3GeV. The values of 
parameters are the same as those in Fig.l. In addition, we considered the 10% asymmetric 
matter fluctuation, a = 0.1. We observe that for E >5GeV, the contribution is much less 
than 1% of the intrinsic T-violation and for lGeV< E <5GeV, the contribution is at most 
3%. Thus we conclude that the asymmetric matter contribution is negligibly small for 
most energies. The small contribution for T-violation asymmetry from the asymmetric 
matter may be understood by observing that the content of the curly parenthesis in 
Eq.(64) vanishes when \kiL\ << 1, \k2L\ << 1 and \kL\ << 1. 

We expect that the linear shape for the asymmetric matter is the biggest deviation 
from the symmetric matter. In order to see the shape dependence of the contribution, 
we consider the cosine shapes and discuss which cosine shape in the Fourier series will 
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give the largest contribution to T- violation. That is, we consider the asymmetric matter 
fluctuation by 

4 

Sa{x) = -anacos{q2n+ix) , (67) 

where g2n+i = (2n + l)7r/L. The coefficient, 4/7r^ is attached by normalizing to the hnear 
shape, {x — L/2)/L = — 4/7r^ X)„=o,i, -(2^ + l)~^cos(2n+ 1)ttx/L. Now we compare the 
contribution for various n by taking a„ = 0.1. We find 



AanaAml^S2xS2yS2zCs f c-^c^-z , SiSz-~z\ . kiL . k2L . kL 

— 7r=;; — — ; ; sm sm sm — 

27r2£;2 \ ki k2 J 2 2 2 

{ k kL cjki kiL s|/c2 k2L] 

L y2n+l ^ '''I y2n+l ^ '^2 y2n+l ^ J 



(68) 



Firstly, we checked that the n = case with ao = 0.1 agrees with the hnear shape 
case with a = 0.1. We confirmed that the difference around edges docs not make any 
difference and the agreement is quite good. Here, we compare n — 1, 2, 3 cases with 
n — case in Fig.4, with the same values of parameters as in Fig.l. That is, we plotted 
{AP)n/ {AP)n=o for n = 1, 2, 3 with = 0.1 for E > 5GeV. We found that as n becomes 
larger, the contribution becomes smaller. The n = 1,2 and 3 cases give about 1/2, 1/5 
and 1/10 times smaller than the n = case, respectively. Thus, we conclude that the 
linear shape case gives the largest contribution to T-violation. 

(b) L = 7332km case 

We use p = 4.21498g/cm^ and the Fourier coefficients, 02, 04, oe and ag are —0.31, 
—0.13, —0.035 and O.Olg/cm^, respectively [9]. In Fig.5, our formula for the sin (5 part 
in Eq.(57) is plotted in comparison with the vacuum contribution. The solid line shows 
the vacuum case, while the 1st and the 2nd order terms are shown by the dotted and 
the dashed fines for 5GeV< E <20GeV, respectively. We observe that the 2nd order 
term is comparable to the 1st order term, and moreover they cancel each other. The 
net contribution shown by the dash-dotted line is quite small. In Fig. 6, the asymmetric 
matter contribution (cos 5 part) is shown in the solid line, in comparison with the sin 5 
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part. As we see from the figure, the cos 5 part is as comparable as the sin 5 due to the 
severe cancellation between the 1st and the 2nd order terms, though it is much smaller 
than the vacuum case. 

Prom the above analysis, we observe the foUowings for L — 7332km: 

(1) We may need to calculate the 3rd order contribution to obtain the accurate formula 
for the sin 6 part in order to check that T- violation is really as small as we obtained. 

(2) The matter fluctuation is no more neglected, because Sa{x)L/2E ~ 1/3 and the 
convergence of the perturbation is not fast. This is because the PREM distribution has 
the simrx/L like structure for L = 7332km and thus the symmetric matter profile can 
not be approximated by the constant (average) density distribution. Therefore, if the 
PREM distribution is correct, the fluctuation from the average medium really gives the 
important contribution for a small quantity such as T-violation. 

8 Summary 

In this paper, we gave the analytical expressions of the contributions from the symmetric 
and asymmetric matter density fluctuations to T-violation in the 2nd order perturbation 
with respect to 5m\^L/2E and da{x)L/2E. Wc found that the contribution to T-violation 
arises from the interference between 6ml^L/2E and da{x)L/2E. The matter fluctuation 
only does not give any contribution to T-violation, which we confirmed by calculating 
up to the 3rd order terms of 5a{x)L/2E. The symmetric and the asymmetric matter 
fluctuations give the effect to the sin 5 and the cos 5 terms, respectively. These analytic 
formula are quite accurate for the distance less than L = 3000km and can be used to 
discuss T-violation analytically. 

By analyzing these formula numerically, we found the following results: Both the 
contributions from the symmetric and asymmetric matter density fluctuations of the 
order of less than 10% to T-violation are small for L = 3000km or a shorter length. 
The use of the average matter density is sufficient for the practical use. Therefore, we 
conclude that the observation of T-violation with L — 3000km or a shorter length will 
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give a quite clear method to determine the CP violation angle, 5. 

For the L = 7332km, the situation changes drastically. We found that the 2nd order 
term from the symmetric matter fluctuation of the order of less than 10% is comparable 
to the 1st order term which includes the effect from the constant (average) matter density, 
and moreover they cancel each other. The net contribution to the sin 5 term becomes 
as small as the contribution from the asymmetric matter fluctuation. This shows that 
the symmetric matter fluctuation becomes important for L = 7332km. This is because 
the PREM distribution has the sinnx/L like structure for L = 7332km and it is not ap- 
proximated by the constant (average) distribution. Therefore, if the PREM distribution 
is correct, the constant (average) density approximation does not work for L — 7332km 
and a longer distance. For L — 7332km, we may need to evaluate the 3rd order term to 
obtain the accurate formula, which is now under the study. 
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Figure 1 : The energy dependence of the T- violation asymmetry for the z/e to channel for 
L = 3000km. (AP)^^ represents T- violation with the symmetric matter profile which is 
proportional to sin 5 as given in Eq.(57), and P = P{i'e ^ ^fi) in oscillation probability in 
the 0-th order. For the symmetric matter profile, the PREM distribution is used. In this 
plot, we use sin 2^^ = sin2^y = 1, sin^^ = 0.1, Am^^ = 5 • lO^^eV^ Am^^ = 3 ■ lO'^eV^ 
and 5 = 7r/4. 
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Figure 2: The energy dependence of T-violation {AP)sg (the dash-dotted hne) in compar- 
ison with the vacuum case (sohd hne) with L = 3000km for (a) 3GeV< E <30GeV and 
(b) lGeV< E <3GeV. We used the PREM distribution for the symmetric matter profile. 
The dashed hne shows the 1st order term, which include the constant (average) matter 
contribution, while the contribution from the 2nd order symmetric matter is indicated 
by the dotted line. The 2nd order term gives negligible contribution and it is hard to see 
in this scale. The difference between {AP)sg and the vacuum case around E =6GeV is 
the matter effect due to the constant (average) density. 
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Figure 3: The energy dependence of the asymmetric matter contribution (a hnear shape) 
to T- violation, which is proportional to cos 5, (AP)c^ with L = 3000km for (a) 3GeV< 
E <30GeV and (b) lGeV< E <3GeV. We assumed the 10% asymmetry of the average 
density. The oscillation parameters are the same as those in Fig.l. By comparing this 
with Fig. 2, we see the asymmetric matter contribution is negligibly small for E > 5GeV 
and about 3% for IGeV <E < 5GeV. 
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Figure 4: The shape dependence of T-violation for the asymmetric matter. Shapes are 
expressed by cos{2n + l)7[x / L . We confirmed that T-violation from the n = case agrees 
with that from the hnear shape case if we choose an appropriate normahzation. We 
plotted the ratios of the contributions from n = 1, 2, 3 and n = 0. We observe that the 
contribution from the higher n is suppressed by about 1/2, 1/5 and 1/10 in comparison 
with the n = case. This shows that the linear shape or = case will give the biggest 
contribution to T-violation. 
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Figure 5: The energy dependence of T- violation from the symmetric matter, (AP)^^ (the 
dash-dotted hne) in comparison with the vacuum case (sohd hne) with L = 7332km 
for 5GeV< E <30GeV. We used the PREM distribution for the symmetric matter. 
The dashed hne shows the 1st order term, which include the constant (average) matter 
contribution, while the contribution from the 2nd order symmetric matter is indicated 
by the dotted line. The oscillation parameters are the same as those in Fig.l. There is 
severe cancellation between the 1st and the 2nd order terms, and the net contribution 
becomes much smaller. This shows that the matter fluctuation gives a sizable effect to 
the sin 6 term and the constant (average) approximation for the symmetric matter does 
not give a good approximation. 
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Figure 6: The comparison between T- violation from the symmetric matter {{AP)sg, 
the dash-dotted hne) and from the asymmetric matter ((AP)c^, the sohd hne) for L = 
7332km. We used the PREM for the symmetric matter and assumed the asymmetry of 
10% in comparison with the average density. The oscillation parameters are the same 
as those in Fig.l. The sin 5 part is not the dominant term in contrast to L = 3000km 
case. This is due to the severe cancellation between the 1st term (including the effect 
from the constant (average) mater) and the 2nd order term from the symmetric matter 
fluctuation. 
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